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Abstract 

If new leptons exist close to the electroweak scale, they can be produced in 
pairs at LHC through standard or new interactions. We study the production of 
heavy lepton pairs in SM extensions with: (i) a Majorana or Dirac lepton triplet, 
as those appearing in type-Ill seesaw; (ii) a lepton isodoublet (N E)l^; (hi) 
a charged isosinglet El^; (iv) a Majorana or Dirac neutrino singlet N and an 

Q ■ additional Z ; gauge boson. It is shown that the trilepton final state which 

has a small SM background, constitutes the golden channel for heavy neutrino 
searches, being very sensitive to Majorana or Dirac neutrinos in triplet, doublet or 
O singlet SU(2)i representations. For higher luminosities, signals in this final state 

can also distinguish lepton triplets from doublets and singlets. The Majorana 
or Dirac nature of the heavy neutrinos is revealed by the presence or not of 
like-sign dilepton £ ± ^ ± signals without significant missing energy Notably, large 

CN ■ £ ± ^ ± signals but with large missing energy are characteristic of Dirac triplets, 

' distinguishing them from the other two models with a heavy Dirac neutrino. 

Further discrimination is achieved with the analysis of the clean £ + l + l~£~ final 

CN ' state. 

in ■ 

o : 
o\ . 

1 Introduction 

> : 

Being a discovery machine, the Large Hadron Collider (LHC) will hopefully uncover 
any new physics close to the electroweak scale. In particular, the several variants of the 
seesaw mechanism proposed to explain light neutrino masses can be tested [1,2]. In 
case that one of these mechanisms is responsible for the light neutrino mass generation 
and seesaw messengers exist around the TeV scale or below, positive signals could be 
observed at LHC. 

There are three types of tree-level seesaw which can originate light neutrino Majo- 
rana masses. The original seesaw [3-6], also known as seesaw I, introduces right-handed 
neutrino singlets N. Seesaw II [7-11] enlarges the SM with a complex scalar triplet A 
with hypercharge Y — 1, which contains three scalars A^, A^, A . Seesaw III [12,13] 
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introduces lepton triplets £ with Y = 0, each containing a charged lepton E and a 
neutral one N. In these three seesaw mechanisms, the lepton number violating (LNV) 
operator [14] 

(05)«,- = tfT&i*4>%L , (i) 

is generated after integration of the heavy degrees of freedom, where L iL are the SM 
lepton doublets with % a flavour index, and is the SM Higgs. Subsequently, when 
electroweak symmetry is spontaneously broken, this operator gives light neutrino Ma- 
jorana masses. The scale at which this operator is generated (i.e. the mass scale of the 
seesaw messengers) is not necessarily very high but it might happen that it is around 
the TeV. In this case, the new leptons predicted in models of type-I and type-Ill seesaw 
might be produced at LHcQ For seesaw III, if the new heavy states E, N have masses 
up to several hundreds of GeV, they can be produced and observed at LHC already at 
the low luminosity phase [1]. For seesaw I, the production of heavy neutrino singlets N 
through SM interactions is very suppressed because of the small heavy neutrino mix- 
ing with SM particles, and signals are unobservable except for relatively "light" masses 
around 150 GeV [21]. However, if new interactions exist, either W' bosons [22-26], 
Z' bosons [27-30] or new scalars [31,32], heavy neutrino singlets can be copiously 
produced, either singly or in pairs depending on the model. 

The observability of the seesaw I III signals and the discrimination among these 
models has been investigated in Ref. [1]. There, an exhaustive analysis of all final 
states was carried out with a complete signal and background calculation, and the 
characteristic features of each seesaw model were highlighted. In this paper we will 
perform a complementary analysis. Our main objective here is to identify the relevant 
signals whose observation or non-observation would discriminate among different mod- 
els with new heavy leptons. In particular, we want to design a strategy to determine 
if new leptons eventually observed at LHC could mediate a type-I or type-Ill seesaw 
mechanism. We will study heavy lepton pair production 

pp — > E + E~ , 
pp -> E ± N , 

pp -> NN , (2) 

where E generically denotes a heavy charged lepton and N a neutral one, in several 
SM extensions with new leptons in different SU(2) £ representations. We consider the 
following additions to the SM particle content: 

1 In seesaw II the scalars predicted can also be produced and observed in this mass range [1] (see 
also Refs. [15-20]). 
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(1) A Majorana lepton triplet £, containing a charged lepton E and a Majorana 
neutrino N. Generically, three such triplets appear in minimal seesaw III reali- 
sations, but as in previous studies [1] we restrict our calculations to the lightest 
one. 

(2) A Dirac lepton triplet E. This is an alternative to the minimal seesaw III in 
which two (quasi-)degenerate Majorana triplets Si, S 2 with opposite CP parities 
form a (quasi-) Dirac triplet £ [2], in analogy with the sometimes called "inverse" 
type-I seesaw with heavy Dirac neutrinos [33-36]. In this case the heavy states 
are two charged leptons E^, E£ and a Dirac neutrino N. 

(3) A lepton isodoublet (N E) LiR , in which the heavy neutrino N is likely to have 
(quasi-)Dirac character. A lepton isodoublet cannot generate the operator in 
Eq. (Q]) but this does not exclude the possibility that they exist, independently 
of the neutrino mass generation mechanism. 

(4) A charged lepton isosinglet E, which can also exist independently of the neutrino 
mass generation. 

(5) A Majorana singlet N and an extra gauge boson Z'. For definiteness we work with 
an E 6 model with a leptophobic Z' x boson and heavy Majorana neutrinos [27], 
assuming that only one of them, with mass < Mz 1 /2, can be produced in 
Z' x decays. But the results obtained here are general because the relative rate of 
the multi-lepton signals produced in pp — > Z' — > NN is only determined by the 
heavy Majorana neutrino decay channels, and the particular type of Z' boson 
considered only affects the total production cross section. Therefore, the results 
shown here can be applied to other models [28, 29] with a trivial rescaling. 

(6) A Dirac singlet N, formed by two Majorana singlets, and an extra Z' boson. We 
work with the same model in Ref. [27] but assuming that two Majorana neutrinos 
form a (quasi-)Dirac neutrino. 

In all the models enumerated above, one or more heavy lepton pairs E + E~, E ± N, 
NN can be produced in hadron collisions. Moreover, in the case of lepton doublets 
and triplets the two states E, N are almost degenerate in mass. But because the pro- 
duction processes and especially the heavy lepton decay channels are different in each 
model, the final state signatures are characteristic and the models can be distinguished 
already at LHC, without the need of precise measurements at a future collider like 
ILC. The model discrimination relies on the simultaneous analysis of different final 
states with two, three and four leptons, and the reconstruction of peaks in invariant 
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mass distributions. We will find that the trilepton signal £^£^1^, with £ = e,/i, is 
common to all models introducing a heavy neutrino§| and has a small SM background. 
Therefore, it constitutes the golden channel for heavy neutrino searches at LHC, in 
particular for the search of seesaw messengers. Moreover, the study of this channel 
can also reveal if a heavy neutrino eventually observed belongs to a triplet or not. 
Its Dirac or Majorana nature can be determined in the like-sign dilepton final state: 
the observation of a £ ± £ ± signal without missing energy clearly indicates its Majorana 
character whereas the absence of such signal points towards one of the other models. 
Like-sign dileptons but with large missing energy can distinguish Dirac triplets, which 
give a large and observable signal, from Dirac doublets and singlets which give much 
smaller ones. Finally, the model identification is completed with the analysis of four 
lepton signals and the search of a resonance in the total invariant mass distribution. 

We must point out that the models considered in this paper do not exhaust all 
possibilities for the addition of heavy leptons and/or new interactions to the SM, but 
rather constitute the cases in which model discrimination seems hardest because the 
heavy leptons are nearly degenerate, they always decay into a light lepton plus a W, Z 
or H boson and their charge cannot be measured in hadronic decays. Two possibilities 
not covered here, in which model discrimination is easier, are: 

• A fourth SM generation of chiral leptons. In this case, the mass splitting be- 
tween the mass eigenstates is expected to be large to be consistent with pre- 
cise electroweak data [37-39], unlike in the cases examined here. Apart from 
this distinctive characteristic, the leading decay of the charged lepton would be 
E — ► W*N, so that E + E~ production would produce events in which the mass 
resonances have more jets than in the models discussed here. Moreover, anomaly 
cancellation requires the presence of new quarks with masses of few hundreds of 
GeV (or new fermions with the same quantum numbers, giving rise to observable 
resonances [39]), which would be produced in pairs and observed already with a 
relatively small luminosity [40,41]. 

• A Majorana neutrino coupling to a new W boson, as for example those appear- 
ing in left-right models. In this case the heavy neutrino can be produced in 
association with a light lepton, 

pp — > W — > £N . (3) 

This process gives signals with up to three charged leptons in some regions of 
parameter space [26], but they can be clearly distinguished from heavy lepton pair 

2 It can also appear in the case of the E singlet but with a very small branching ratio. 
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production by the reconstruction of only one heavy resonance which, together 
with an energetic light lepton, produces a Jacobian peak at the W mass. 

In this work we do not consider scalar triplet production (seesaw II) either, which gives 
signals kinematically very different with sharp peaks in like-sign dilepton invariant 
mass distributions. It is interesting, however, to note that in seesaw II the trilepton 
signals are the most important as well [1]. Multi-lepton signals can also appear in pair 
production of new quarks [42] but they can be easily distinguished from the models 
studied here by mass reconstruction and also by the presence of significant single lepton 
signals with several fe-tagged jets. 

The structure of this paper is as follows. After this introduction, we briefly sum- 
marise in section [2] the relevant Lagrangian terms, the production processes at LHC 
and the decay channels of the heavy leptons in the models studied. In section [3] we 
describe how the signals and backgrounds are simulated. The results for the £ ± £ ± £ T , 
and final states are presented in sections 0HH1 Section 7 summarises 

these results and shows explicitly how the six models considered in this work can be 
distinguished. We draw our conclusions in section [H In appendix lAl we give the partial 
widths for heavy lepton decays in the different models studied. 



2 Description of the models 

Here we present the different models studied in turn, collecting the interactions relevant 
for our analysis and enumerating the production processes at LHC and the allowed 
decay channels of the heavy leptons. At the end of this section we summarise the main 
features and compare among the different models. 

2.1 A Majorana lepton triplet 

The relevant Lagrangian for any number of Majorana triplets Sj has been previously 
given in Ref. [1], and we follow the notation in that work. For a single triplet S, the 
interactions of the heavy leptons E, N with SM leptons I, v\ is, at first order in the 
light-heavy mixing Vin, 

C w = -g (E-fN W~ + N^E W+) 

- ±* (v lN I^PlN w~ + v; N nyplI w+) 

- g (V lN E^P m W~ + V t * N Vi^PrE W+) , 
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C z = gc w E 1 »EZ i 



ZCyy 

+ -ji— (v lN h»p L E + v; N e^PlI) z, , 

C 1 = eE-fEA ll , 

+ {Vin WrE + V; N EP L l) H . (4) 

Heavy lepton pairs are produced by the gauge ZEE, jEE and WEN interactions, 

qq -> Z*h* -> E + E~ , 

qq' -> -> £± JV , (5) 

(note that there are no ZNN interactions because the triplet has zero hypercharge, 
and iViV pairs are not produced). The heavy leptons E, N decay to SM leptons plus 
a gauge or Higgs boson: 

E+^vW + , E + ^l + Z, E + ^l + H, 

N^TW + , N^l + W~, N -> vZ , N -> i/# . (6) 

The two heavy states are nearly degenerate in mass, with a small splitting due to 
radiative corrections, and the decays between the heavy states are very suppressed [43]. 



2.2 A Dirac lepton triplet 

A (quasi-)Dirac lepton triplet is formed by two (quasi-)degenerate Majorana ones with 
opposite CP parities. As it was shown in Ref. [2], the heavy fields can be redefined in 
such a way that the Lagrangian is written in terms of two charged leptons E± , E% (the 
fermion is positively charged in the second case) and a Dirac neutrino N, and lepton 
number is conserved up to effects of the order of light neutrino masses. At first order 
in Vim we have 

C w = -g( E^N - N^Et) W; - g (N-fE^ - E+>fN) W+ 
- JL (y lN I^PlN W- + V; N N^PlI w+) 

+ g {Vin vii^PLEi w~ + v; N Et-fPm w;) , 
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c z = 9 c w ( ei^ei - E+>fE+) 

+ ^ {Vin viYPlN + v* N N^Pm) z, 

+ -JL- (v lN i-fp L Ei + v; N e^PlI) z, , 

£ 7 = e ( E^-fEZ - Et^Et) A, , 

+ {Vin Tp r ei + v lN e- x p l i) H . (7) 

Heavy lepton pairs are produced by the gauge ZEiEi, r yE i Ei and WEiN interactions, 

qq -> Z*/j* - E+Er , 

qq'^W*^E±N, (8) 

Since there are two charged fermions instead of only one (and a Dirac neutrino is 
equivalent to two Majorana ones), the total heavy lepton production cross section 
is twice larger than for a Majorana triplet. The decays have some differences with 
respect to a Majorana triplet because E{ does not couple to light neutrinos, does 
not couple to light charged leptons and N decays conserve lepton number. Thus, the 
allowed ones are 

El -> l~Z , El -> l~H , E+ -> vW + , 

N^TW + , N -> vZ , N -> uH . (9) 
2.3 A lepton isodoublet 

The Lagrangian for a lepton isodoublet can be found in Ref. [44]. With our notation, 
the terms involved in heavy lepton production and decay are, at first order in Vin, 

9 



C w = -^=( E-fN W~ + N-fE W+) 

- JL (y lN T-fP R N W- + V lN N^PrI w+) 



C z = ([-1 + 2s$y]E>fE + N-fN) Z, 

+ 7T- ( v i* WPrE + V lN E^PrI) Z, , 

ZC\y 

£ " = ( V ™ IP l E + V i*n EP R l ) H • ( 10 ) 
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Note that the neutrino is likely to be a Dirac fermion because the renormalisable 
gauge-invariant doublet mass term 

£ mass = ~ m D L 4R L 4L + H.c. , (11) 

with L 4 = (N E) T , implies a Dirac mass. Additional Majorana masses may appear 
from dimension-five operators, 

£Lm = x ^ftu L + LfoirpLut . (12) 

However, if the physics generating these operators is the same as the one yielding the 
light neutrino mass operator in Eq. ([]]), one would expect that the Majorana mass 
terms, of order CijV 2 /A, are much smaller than mr>. In this case, N is a (quasi-)Dirac 
fermion. 

The production processes are the same as for a Majorana triplet but now neutral 
pairs are also produced because they couple to the Z boson, 

qq -> Z*/ 7 * -> E + E- , 
qq -> W* -> E ± N , 

qq^ Z* -> NN . (13) 

A further difference is in the size of the couplings, e.g the coupling to the W boson is 
reduced by a factor a/2, and so the E ± N cross section is a factor of two smaller. In 
this model the heavy lepton decays are also different, being allowed only 

E + ^l + Z, E + ^l + H, 

N^TW + . (14) 



2.4 A charged singlet 

The Lagrangian for a charged lepton isosinglet can also be found in Ref. [44]. In our 
notation, and at first order in the light-heavy mixing V Evi , the relevant terms are 

C w = (V Evi EYPlVi W- + V Evi vrfP L E W+) , 

£z = - 9 -^E^EZ, 

c w 

+ J- (V Eui E-fPrf + V Evi T-fP L E) Z, , 
C y = eE'fE Ap, 

Ch = ^ EPl1 + VL < IPrE "> H (15) 



w 
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In this model only charged lepton pairs are produced, 



qq -> Z*h* -> E + E~ , (16) 
which later decay in the three possible modes 

E+^vW + , E + ^l + Z E + ^l + H. (17) 

2.5 A Majorana neutrino and a Z' boson 

A heavy Majorana neutrino which is a singlet under SU(3) x SU(2)l x U(l)y interacts 
with SM fields via a small mixing Vjjv with SM fermions (for a detailed derivation of the 
Lagrangian see for example Ref. [45]). Its interactions are therefore suppressed, being 
at least of order Vw (ZNN interactions are of order V t 2 N ). However, the interactions 
with an extra Z' boson are not suppressed if N is not a singlet under the Z' gauge 
group U(l)', and are determined by the heavy neutrino charge Q under this extra U(l)'. 
At first order in Vw, the relevant Lagrangian is 

Cw = {V lN T-fP L N W~ + V{ N NYPlI W+) , 

c z = ~^-^{v in p l -vc n p r )n z„, 



2c 



w 



C-h = -^udVwPn + V^NH, 

C z , = -g'^N^NZ'^. (18) 

In the model considered here [27] we have Q = — 3, and we have generically denoted 
the U(l)' coupling as g'. We assume that it equals the SM coupling gy of U(l)y, but 
this may change by renormalisation group evolution effects. Heavy neutrino pairs can 
be produced with the exchange of an s-channel Z' boson, 

qq -> Z' -> NN , (19) 

and the heavy neutrinos decay giving a W, Z or H boson plus a light lepton, 

N^TW + , N^l + W~, N^vZ, N^vH. (20) 



2.6 A Dirac neutrino and a Z' boson 

Finally, in the case of a (quasi-)Dirac neutrino composed by two (quasi-)degenerate 
Majorana neutrinos with opposite CP parities and equal charges Q under U(l)', the 
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relevant Lagrangian is 




(21) 



Heavy Dirac neutrino pairs can be produced in the same process 



qq 



NN , 



(22) 



and the cross section is twice larger than for a Majorana N because the symmetry 
factor for identical particles is not present in this case (the Z' — ► NN width is also a 
factor of two larger) Jf| Heavy neutrino decays are the same except for the absence of 
the LNV one, 



2.7 Summary 

In each of the preceding subsections we have indicated the heavy lepton pair production 
processes present in each model. For better comparison of the different models, we 
summarise in Table Q] the coupling constants appearing in the production vertices 
involved, coupling two heavy leptons to a gauge boson. From now on we will use the 
labels Em, E d to refer to the models with Majorana and Dirac triplets, respectively, 
NE d and E s for the doublet (NE)l,r and singlet El,r, and Z'N M , Z'N^ for the 
SM extensions with a leptophobic Z' x boson and a Majorana or Dirac neutrino. The 
production cross sections are presented in Fig. [Q as a function of the heavy lepton mass 
(as well as the Z' mass in the corresponding models). The decay of the heavy leptons 
takes place in the channels indicated, with partial widths collected in appendix [A] 
for reference. Nevertheless, the important quantities for LHC phenomenology are the 
relative branching ratios. Summing over light leptons I = e, //, r, the branching ratios 
into W, Z and H bosons are independent of the mixing. Their values for m E>N ^> 
Mw, Mz, Mh, are collected in Table El This table illustrates the important differences 
among the models considered, which make the discrimination based on multi-lepton 

3 Note that the Feynman rule for Majorana fermions contains an extra factor of two to account for 
the two possible Wick contractions, so that the Z'NN vertex is —ig'Q^^^h as in the Dirac case. 
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(23) 



10 





WEN 


ZEE 


jEE 


ZNN 


Z'NN 




9 


gc w 


e 





— 




9 


gc w 


e 





— 


E D , E£ 


9 


gc w 


e 







NE d 


9/V2 


g/c w (s 2 w - 1/2) 


e 


9/ (2cw) 




E s 




gsw/cw 


e 






Z'N M 








~ 


g'Q/2 


Z'N B 








~ 


g'Q 



Table 1: Coupling constants of the gauge vertices involved in heavy lepton pair pro- 
duction, for the different models considered (the notation is given in the text). The 
labels E±, E^ refer to the two charged leptons in the model with a Dirac triplet. A 
dash indicates that some of the particles in a vertex does not exist in a given model. 
A zero indicates that the particles exist but they do not couple. 
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Figure 1: Left: Cross section for heavy lepton pair production in different models. 
For the Dirac triplet (S D ) the cross sections are twice larger than for the Majorana 
one (S M ). Right: cross section for heavy Majorana neutrino pair production from Z' x 
decays. For Dirac neutrinos the cross section is two times larger. 



signals quite effective. Note that for masses m EiN = 300 GeV as considered in the 
following sections, the branching ratios are already close to the values presented in this 
table. 
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Table 2: Decay branching ratios of the heavy leptons in the limit ttle^n 3> 
M w , M z , Mj{, for the different models considered (the notation is given in the text). 
The labels E± , E% refer to the two charged leptons in the model with a Dirac triplet. 

3 Multi-lepton signal generation 

The analysis pursued in this work, aiming to discriminate among several models all 
giving various multi-lepton signals in different decay channels, is somewhat demanding 
from the point of view of the simulation. It requires to generate all the signal con- 
tributions because many different heavy lepton decay channels, with the subsequent 
W/Z/H decay, can lead to the same charged lepton multiplicities. The complete sig- 
nal generation has been done with the program Triada [1] extended to include the 
models with an isodoublet (N E)l,r, an isosinglet El,r and Z' production with decay 
to a Majorana or Dirac neutrino. All the signal processes enumerated in the previous 
section, with all the possible decays of the heavy leptons E, N and the W/Z bosons, 
are included. The Higgs boson decay, which does not carry any spin information, is 
left to the parton shower Monte Carlo. Signals have been generated with statistics of 
300 fb -1 and rescaled to a reference luminosity of 30 fb _1 , in order to reduce statistical 
fluctuations. The SM background, consisting of the processes in Table [3l is generated 
using Alpgen [46] with a Monte Carlo statistics of 30 fb -1 . Additional SM processes 
which were previously shown to be negligible after selection cuts [1] are ignored in this 
work. Signals and backgrounds are passed through the parton shower Monte Carlo 
Pythia 6.4 [47] to add initial and final state radiation (ISR, FSR) and pile-up, and 
perform hadronisation. For the backgrounds we use the MLM prescription [48] for the 
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matching to avoid double counting between the matrix-level generator and the parton 
shower Monte Carlo. We use the fast simulation AcerDET [49] which is a generic LHC 
detector simulation, neither of ATLAS nor of CMS, finding good agreement with our 
previous results [1]. 



Process Decay Events 

tinj, n = 0, . . . , 5 semileptonic 6.1 M 

tinj, n = 0, . . . , 5 dileptonic 1.5 M 

tW all 1.6 M 

Wtinj, n = 0,...,4 W -► Iv 5.1 K 

Zbbnj, n = 0, .... 4 Z -► Z+r 200 K 

Zttnj,n = 0,...,4 Z Z+Z - 1.87 K 

Wnj, n = 0, ...,3 W-»Zi/ 290 K 

WZry, n = 0,...,3 W Zi/, Z -► Z+Z~ 37.7 K 

ZZnj, n = 0, . . . , 3 Z -► Z+Z~ 3.74 K 

WWWnj, n = 0,...,3 2W Zi/ 1.47 K 



Table 3: Background processes considered in the simulations, with nj standing for n 
additional jets at the partonic level. The second column indicates the decay modes 
included (where / = e, //, r). The last column corresponds to the number of events 
after matching for a luminosity of 30 fb _1 , with K and M standing for 10 3 and 10 6 
events, respectively. 

We remark that the use of at least a fast simulation of the detector is essential for 
this and other similar studies, because some of the most important SM backgrounds 
have charged leptons resulting from b quark decays, and these cannot be estimated in 
a parton-level analysis. In fact, the most recent analyses for supersymmetry searches 
performed with a full detector simulation [50] confirm the well-known fact that tinj 
is probably the largest (and most dangerous) SM source of like-sign dileptons. From 
the point of view of the signal, the use of a fast detector simulation is also necessary 
because some of the contributions to the multi-leptonic final states studied arise when 
more charged leptons are produced but missed by the detector. 

In the following sections we will take heavy lepton masses m EjN = 300 GeV for our 
simulations. Production cross sections are independent of the mixing and, except for 
unnaturally small values, the heavy leptons will decay well inside the detector. Note 
that the heavy lepton widths are larger than the b quark width for V > 10 -7 , and 
present indirect constraints are of order V < 10 _1 [51]. In definite seesaw models, 
heavy lepton mixing is also related to light neutrino masses. However, we do not 
assume any particular model-dependent relation between light neutrino masses and 
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heavy lepton mixing. Instead, we assume that heavy leptons only couple to the first 
generation, bearing in mind that for an arbitrary mixing with the first two generations 
the results obtained would be approximately the same (the signals are equivalent and 
at high transverse momenta the SM backgrounds involving electrons and muons have 
roughly the same size). On the other hand, if E, N only mix with the tau the signals 
would be very difficult to observe [1]. For the Z' x boson we will conservatively take 
a mass of 650 GeV, because this is approximately the location of the maximum in 
the heavy lepton pair invariant mass distribution for m,E,N = 300 GeV when they are 
produced by off-shell W/Z bosons. Hence, for such mass the identification of an s- 
channel resonance would be more difficult. This Z' x mass is not excluded by present 
Tevatron searches for it resonances. For the Z' x model considered, the cross section 
into tt final states (assuming g' = gy) is about 0.2 pb, well below the 95% confidence 
level limit for this mass, which is around 0.6 pb [52]. The Higgs boson mass is taken 
as M H = 115 GeV, in which case it mainly decays into two jets, H — > bb,cc,gg and 
seldom produces leptons, only when H — > t + t~ with r leptonic decay. For a heavier 
Higgs decaying into W + W~ , ZZ , the multi-lepton signals examined here would still be 
present but some signals with higher lepton multiplicity, originating from heavy lepton 
decays involving a Higgs boson and H — > W + W~ , ZZ, would also be present and might 
be of interest. In addition, several SM backgrounds (for example, W ± H — > W ± W + W~ 
plus jets) would be enhanced. A dedicated analysis is required to examine the precise 
discovery potential of each channel in such case. 

Finally, it is worth commenting here about the statistical prescriptions used to 
determine a possible discovery. If the background can be precisely known or directly 
estimated from data, for instance if the signal shows up as a sharp peak, the statistical 
significance is Sq = S/ \fB, where S, B are the number of signal and background events, 
respectively. For small B, this estimator is replaced by the P-number using Poisson 
statistics. The discovery criteria used in this work are: (i) statistical significance larger 
than 5<r, and (ii) the presence of at least 10 signal events. In most of our results 
the luminosity required for discovery does not depend on a very precise background 
normalisation because the background is tiny, so that the signal significance is well 
above 5a and the limit is determined by having 10 signal events. 
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4 Final state £±£±1* 

This final state is the most characteristic one of a heavy neutrino, and it can appear 
both in E^N and NN production, for example in the decay channels 

E+N ^l+Zi+W*, Z -> qqjvv.W ' -> Iv , 

E + N ^ £ + H£ ± W T , H^qq,W^£v, 

nn -> , ww -> ^ , (24) 

irrespectively of the Dirac or Majorana character of N. For a lepton triplet, it can also 
be produced in several other decays, for example 

E + N -> t + ZvZ, ZZ ^ l + T qq/vv , 

E + N -> tZuH/tHuZ , Z ^ £ + r,H ^ qq. (25) 

Note that iV — >■ z/Z, iV — > z/if do not take place in the case of a lepton doublet, 
which can be exploited to distinguish them from triplets. With this aim, we perform 
two complementary analyses for the trirlfi final state, dividing the sample into two 
disjoint ones. First, we perform a generic analysis as in Ref. [1], which can detect 
the presence of a heavy neutrino singlet, doublet or triplet decaying in the channels of 
Eqs. (J24j). In this analysis we reject events with a Z boson candidate, that is, with an 
opposite-charge lepton pair with an invariant mass consistent with Mz- This sample is 
labelled as £+£+£+ (no Z). In second place, we perform a new specific analysis to search 
for the decays in Eqs. (1251 and determine whether the neutrino belongs to a triplet. 
In this case we only accept events with a Z boson candidate. This sample is labelled 
as £+£+£+ (Z). The common event pre-selection criteria for the two analyses are: (i) 
the presence of two like-sign leptons l\ and £2 with transverse momentum pt > 30 
GeV and an additional (and only one) lepton of opposite sign, with px > 10 GeV; (ii) 
two hard jets with px > 20 GeV. The cut on transverse momenta of the like-sign pair 
greatly reduces the SM background from tinj production. The requirement of two 
jets is used in the kinematical reconstruction, and also reduces the WZnj background. 
We must point out that this background, as simulated with Alpgen, does not include 
off-shell photons and uses the narrow width approximation for both bosons. The effect 
of off-shell photons is important in general [53] but in this case this contribution is 
reduced by the high-p^ requirement on charged leptons. More important is the effect 
of the Z width: the narrow approximation underestimates the background for the 'no 
Z' sample and overestimates it in the complementary l Z' one. With a comparison of 
the reconstructed £ + £~ distributions for Z/j* and Z on-shell after detector simulation, 
we estimate that for the 'no Z 7 trilepton final state the background can be at most two 
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times larger than the values given here. This can be compensated, at the cost of some 
signal efficiency loss, by a wider interval for rejection of events with a Z candidate, and 
the discovery potential would be very similar. Model discrimination, of course, would 
not be affected. 

4.1 Final state PlH* (no Z) 

In this analysis we ask for event selection the absence of a Z boson candidate: neither 
of the two opposite-sign lepton pairs can have an invariant mass closer to Mz than 10 
GeV. This is very useful to remove WZnj production but also eliminates several of the 
signal channels. We collect in Table H] the number of signal and background events for 
each process and model after pre-selection and selection cuts. 
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Table 4: Number of events in the t^i^t^ (no Z) sample for the signals and main 
backgrounds with a luminosity of 30 fb -1 . 

The event reconstruction is performed in three steps, following this procedure [1]: 

1. The momentum of the Z or H boson decaying hadronically is reconstructed as 
the sum of the momenta of the leading and sub-leading jets. 

2. One of the heavy charged leptons L (E or N, depending on the process) can 
be reconstructed from this boson and one of the two like-sign leptons, and the 
heavy neutrino N from the two remaining leptons (with opposite charge) and 
the missing neutrino momentum. The longitudinal component of the neutrino 
momentum is neglected for the moment, and the transverse component is taken 
as the missing energy. There are two possibilities for this pairing, and we choose 
the one giving closest invariant masses for the reconstructed L and N. 
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Figure 2: Left: reconstructed heavy lepton masses for the signals in the t^l^fO* (no 
Z) sample. Right: the same for the SM background and the background plus the 
Majorana triplet signal. The luminosity is 30 fb _1 . 



3. The N reconstruction can be refined by including the longitudinal neutrino mo- 
mentum. We select among the two charged leptons the least energetic one £ s , 
and require that its invariant mass with the neutrino is Mw, 

( P i s +Pu) 2 = M^, (26) 

taking the transverse components of p u as the missing energy. This quadratic 
equation determines the longitudinal neutrino momentum up to a twofold am- 
biguity, which is resolved selecting the solution with smaller (p v ) z . In case that 
no real solution exists, the transverse neutrino momentum used in Eq. (1261) is 
decreased until a real solution is found. 

The reconstructed heavy lepton masses after selection criteria are shown in Fig. [2] 
(left) for the five models which give an observable signal. We point out that one of the 
resonances, labelled as N, can be identified as being a neutral lepton, because the peak 
appears in the invariant mass distribution of two opposite-charge leptons plus missing 
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energy. The identity of the other resonance, labelled as L, cannot be established be- 
cause the charge of the hadronic jets cannot be measured. In fact, depending on the 
process and the model, it can be a charged or neutral heavy lepton (see Table [4]). In 
Fig. [2] (right) we show the same distributions for the SM background and the back- 
ground plus the Majorana triplet signal, which is the smallest one. We define the peaks 
as the intervals 

240 < m"l c < 360 GeV , 

240 < m^ c < 360 GeV , (27) 

and show in Table H] the number of signal and background events after these cuts. The 
statistical significance for 30 fb _1 of the relevant signals (neglecting the background 
uncertainty) is collected in Table [5j with the luminosity required to achieve 5cr discov- 
ery. Due to the smallness of the background, the discovery luminosity is determined 
in each case by the requirement of having at least 10 signal events. For such small 
luminosities the SM background is tiny, around one event or less, and a very precise 
normalisation is not very important. Notice also that the signals (and thus the sta- 
tistical significance) are larger for a lepton doublet than for the triplets, despite the 
smaller production cross section. This is mainly due to the larger (100%) branching 
ratio of the mode N — ► £~W + , and also to the presence of NN production. 





So 


L 




So 


L 




24.6 


3.3 fb" 1 


Z'N M 


38.3 


2.1 fb" 1 




54.8 


1.5 fb" 1 


Z'N B 


76.3 


1.1 fb- 1 


NE A 


74.0 


1.1 fb" 1 









Table 5: Statistical significance of the relevant signals for 30 fb 1 , and luminosity L 
required to have a 5cr discovery in the E^i^P (no Z) sample. 

With the results shown, we observe that this final state is very sensitive to the 
presence of heavy neutrinos in all the models considered (except the one with the 
charged singlet, which does not have a neutral lepton) Q We can then ask ourselves 
whether one could already distinguish some of these models. This is indeed possible, 
although model discrimination probably requires more luminosity than 5cr discovery. 
By examination of the heavy lepton pair invariant mass, one can determine if these 
leptons are produced by the exchange of an s-channel resonance, as it is the case in 



4 Note that the sensitivity to Majorana triplets can be improved with a more inclusive analysis 
without event reconstruction, and the luminosity required for 5a discovery can be reduced to 1.7 



fb- 1 [1]. 
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Figure 3: Left: invariant mass distribution of the heavy lepton pair in the i^i^P 
(no Z) final state, after reconstruction cuts, normalised to 30 fb -1 . Right: possible 
experimental result for 10 fb -1 , where the numbers of signal events are 31 (Em), 68 
(S D ), 92 (NE d ), 48 (Z'N M ) and 95 (Z'N D ). 



the models with an extra Z' x boson. We show in Fig. [3] this kinematical distribution 
for the signals only, separating for clarity the models in which the heavy neutrino has 
Majorana (up) and Dirac nature (down), which can be distinguished by other means 
(see section 15 . lh . On the right we show possible experimental results for 10 fb -1 , 
obtained by making random fluctuations with a Poisson distribution of the bins in the 
distributions, and normalising to the total expected number of events. The background 
is much smaller than the signals and has not been included. It is quite clear that the 
presence of a resonance can be detected or excluded, possibly with a smaller luminosity, 
but we do not address this issue quantitatively. 

4.2 Final state tl^F (Z) 

The trilepton sample with a Z boson candidate suffers from a large WZnj background. 
Nevertheless, several of the decay channels in Eqs. (1251) produce a sharp peak in the 
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Table 6: Number of events in the (Z) sample for the signals and main back- 

grounds with a luminosity of 30 fb -1 . 
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Figure 4: Left: reconstructed heavy charged lepton mass for the signals in the £ £ £ T 
(Z) sample. Right: the same for the SM background and the background plus the 
Majorana triplet signal. The luminosity is 30 fb -1 . 



trilepton invariant mass distribution, corresponding to the heavy charged lepton mass 
m E . This peak can be observed over the large background. As event selection criteria 
we require: (i) a Z boson candidate, with two opposite-charge leptons having an invari- 
ant mass between Mz — 10 GeV and Mz + 10 GeV; (ii) missing energy j/t > 30 GeV. The 
number of signal and background events after event pre-selection and selection is given 
in Table [61 The heavy E mass is simply reconstructed as the three-lepton invariant 
mass. Its distribution is shown in Fig. [4] (left) for the five non-negligible signals. For 
the lepton triplet the peaks could be seen over the SM background, as it is shown in 
the right side of this figure for the Majorana triplet. These results can be understood 
with the following considerations: 
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So L Sq L 

S M 5.5 25 fb" 1 S D 6.7 16.6 fb" 1 

Table 7: Statistical significance for 30 fb -1 of the relevant signals, and luminosity L 
required to have a 5<r discovery in the i^i^P (Z) sample. 

1. The off-peak contributions in the lepton triplet models result, for example, from 
the decay channels in Eqs. (1251) with N — > vZ — > v£ + £~. In this case the 
resonance would be seen in the invariant mass of two jets plus the other lepton, 
all of them being E decay products. However, the corresponding peak is very 
broad due to the worse jet energy resolution, and difficult to see over the large 
WZnj background. 

2. As expected, in processes where a charged lepton is not present, as in Z' x — > NN, 
the distributions do not display a peak. 

3. For a lepton doublet the decay channel EN — > £Z£'W, with £' missed by the 
detector and Z — > £ + £~, W — > qq, would also produce a trilepton signal with a 
peak at m^. However, this process is removed by the missing energy requirement, 
and the rest of the decay channels do not have any resonance in the trilepton 
invariant mass, as it can be seen in Fig. [H This feature clearly distinguishes a 
lepton doublet from a (Majorana or Dirac) triplet. 

Defining the peak region as 

280 < m1 c < 320 GeV , (28) 

and performing a kinematical cut on this reconstructed mass, we obtain for each process 
the number of events listed in Table [61 The statistical significance of the signals and 
the luminosity required for 5a discovery are presented in Table We assume that the 
background can be accurately normalised with off-peak measurements of the trilepton 
invariant mass distribution, as it is apparent in Fig. 0] (right). It is clear that this 
sample improves little the global statistical significance of the trilepton channeljfl But 
the observation of a positive signal here, with about half the size of the signal in the 
sample without Z candidates, is a proof of the triplet nature of the heavy neutrino. And 

5 The sensitivity shown here may be significantly improved by further reducing the background 
if we ask for the presence of two hard jets and set a more stringent cut on m Y ^ c . However, in a 
real experiment it may be difficult to spot the presence of a peak with small statistics due to the 
background fluctuations, and to be conservative we have not performed these optimisations. 
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it also proves the presence of a charged heavy lepton, although this can be established 
more easily in the four lepton final state, as it will be shown in section El 

5 Final state tt^ 

The like-sign dilepton final state without significant missing energy is characteristic of 
a heavy Majorana neutrino whose decays violate lepton number. It can be produced 
in decays such as 

E + N -> £ + Z £ + W~ , Z -»• qq, W -»• qq' , 

E + N -> £ + H £ + W~ , H ^qq,W -> qq' , 

NN -»■ £±W T £±W T , WW -> gg'gg' . (29) 

Large like-sign dilepton signals can also arise from lepton number conserving (LNC) 
processes, most notoriously in 

E+N -> vW+ £+W~ , W + -> £+v, W~ -> gg (30) 

in the Dirac triplet model, being lepton number balanced by final state neutrinos 
(remember that E + is a fermion, not an anti-fermion). They can also appear when 
more charged leptons are produced but are missed by the detector, for example 

E + N -+£+Z£-W + , Z -> qq,W -> 

NN^£ + W~rW + , W~ qq,W + ^ £ + v , (31) 

with t~ missed by the detector. The presence of such signals should not constitute 
a surprise, since several SM processes, for instance W+W+nj and WZnj production, 
give like-sign dileptons in this way. Hence, as in the trilepton channel we will perform 
two different analyses with disjoint event samples. In the first one we will, as in 
Ref. [1], impose selection criteria to isolate the truly LNV processes characteristic of 
a Majorana fermion. In particular, we will require the absence of significant missing 
energy. In second place we will perform a new analysis for events with large missing 
energy, aiming to isolate the large signal from the decay channel in Eq. fl30l . The 
common pre-selection will be to ask the presence of two like-sign leptons £\ and £2 with 
transverse momentum pt > 30 GeV. 

5.1 Final state e iz £ ± (no jfr) 

In order to keep only the LNV signals characteristic of a Majorana fermion, we ask for 
event selection (i) missing energy pr < 30 GeV, and (ii) the presence of at least four 
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Table 8: Number of events in the t^t^ (no p/r ) sample for the signals and main 
backgrounds with a luminosity of 30 fb -1 . 

hard jets with pt > 20 GeV. In Table [8] we collect the number of signal and background 
events at pre-selection and selection. We observe that our selection cuts very efficiently 
select only the signals with LNV neutrino decays, rejecting for example a 97.5% of the 
large Dirac triplet signal. The event reconstruction is performed as follows [1]: 

1. We associate each charged lepton to a pair of jets in all possible ways, using the 
four jets with larger pt- 

2. Among the six possibilities, we choose the one minimising the difference between 
the two jj and the two £jj invariant masses, 

(m jlj2 - m M J 2 + (m £ljlj2 - m l2jzji f . (32) 

Note that for the leading signal contributions two of the jets in principle corre- 
spond to a hadronic W decay and the other two to a Z or Higgs boson decay. 
However, if a wrong assignment is made, it is expected that the invariant mass 
differences will be larger. 

We present in Fig. [5] (left) the reconstructed heavy lepton masses for the two models 
with heavy Majorana neutrinos which are the only ones yielding observable signals. 
On the right side we show the same distributions for the SM background and the 
background plus the Majorana triplet signal, which has a similar size as the Majorana 
neutrino singlet one. The presence of new resonances is apparent and the background 
normalisation should not be a problem using data away from the peak regions. We 
have labelled the resonances as Li, L 2 because their identity cannot be established 
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Figure 5: Left: reconstructed heavy lepton masses for the signals in the £ ± £ ± (no 
j/r ) sample. Right: the same for the SM background and the background plus the 
Majorana triplet signal. The luminosity is 30 fb -1 . 



(their charge is not measured). Actually, for the Majorana triplet signal they are in 
most cases a charged lepton and a neutral one, while for the Z' x model both resonances 
are neutral (see Tabled]). Therefore, the observation of this signal alone points towards 
the existence of a Majorana (and hence electrically neutral) fermion but its charge is 
not directly measured as in the trilepton final state. 

We define the peak regions as 

250 < m£ c < 350 GeV , 

250 < m'H < 350 GeV , (33) 

and show in Table [8] the number of signal and background events after these cuts. The 
statistical significance of the signals for 30 fb -1 (neglecting the background uncertainty) 
is collected in Table [U with the luminosity required to achieve 5<r discovery. Due to the 
smallness of the background, the discovery criterion is again controlled by the require- 
ment of having at least 10 signal events, and an accurate background normalisation is 
not crucial because for the discovery luminosities it is very small, around one event. 
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So L 
S M 35.0 2.1 fb" 1 



So L 
Z'N M 38.3 2.3 fb" 1 



Table 9: Statistical significance for 30 fb 1 of the relevant signals, and luminosity L 
required to have a 5cr discovery in the (^(^ (no fa) sample. 



We also address here the discrimination of the two models with M ajorana neutrinos 
in the like-sign dilepton final state. This can easily be done with the heavy lepton pair 
invariant mass distribution, shown in Fig. El which displays a peak for the production 
by an s-channel Z' boson and is more flat for the lepton triplet. In the same figure we 
show on the right panel a possible experimental result for 10 fb -1 , obtained with random 
fluctuations with a Poisson distribution of the bins of the distributions, normalised to 
the expected total number of events. Both models give quite different results, and the 
discrimination seems easy once that sufficient luminosity is collected. 
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Figure 6: Left: invariant mass distribution of the heavy lepton pair in the 
j/t) sample, after reconstruction cuts. Right: possible experimental result for 10 fb -1 . 



5.2 Final state (j/ T ) 

The selection criteria to isolate the LNC dilepton signals from Dirac triplets are: (i) 
large missing energy p'x > 50 GeV; (ii) a transverse momentum > 100 GeV for the 
leading charged lepton; (iii) at least two hard jets with pt > 20 GeV. The number 
of signal and background events is gathered in Table [101 Notice that the Majorana 
triplet also produces a large signal, from a decay similar to the one in Eq. (1301) but with 
lepton number violation (in this model E + is an anti-fermion). The triplet signals are 
reconstructed as follows: 
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Table 10: Number of events in the £ (pr ) sample for the signals and main back- 
grounds with a luminosity of 30 fb _1 . 

1. We first reconstruct the W mass as the invariant mass of the two jets with largest 
transverse momentum. The result is shown in Fig. [7] (up). 

2. There are two possible pairings with the two charged leptons, iijj and ^33, to 
reconstruct the heavy lepton invariant mass. We use both to construct a plot with 
two entries per event shown in the same figure (middle), which clearly displays 
a peak at the true mass for the triplet signals, and small peaks also for the rest. 
This plot can be used to determine the resonance mass. 

3. The reconstruction is completed selecting among the two possibilities the one 
giving a reconstructed mass close to the value determined from the peak in the 
previous distributions. This choice introduces some bias in the non-triplet signals 
and the SM background, as seen in Fig. [7] (down) but keeps the signal larger. 

Notice that despite the appearance of small peaks in the non-triplet signals the dif- 
ferences with the lepton triplet are huge. We perform kinematical cuts around the 
peaks, 

60 < m^ c < 100 GeV , 
260 < m r l c < 340 GeV , (34) 

which are sufficient to reduce the background. The corresponding numbers of events are 
given in Table UM The statistical significance of the signals, ignoring the background 
uncertainty, is given in Table El with the luminosity required to have a 5a significance. 
For the Majorana and Dirac triplet signals, which are quite large and observable already 
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Figure 7: Up: reconstructed W mass for the signals (left) and the Dirac triplet signal 
with background (right) in the £ ± £ ± {p/p ) sample. Middle: the same for the £ijj 
and i-ijj invariant mass distribution (two entries per event). Down: the same for the 
reconstructed heavy lepton mass. The luminosity is 30 fb _1 . 
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with very small luminosity, the background normalisation is not crucial. In the rest 
of models the figures in Table [TT1 are optimistic but this even enforces our argument, 
which is to point out that lepton triplet signals can be seen in this final state while for 
the rest of models the signals are much harder to observe. 
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Table 11: Statistical significance for 30 fb 1 of the relevant signals, and luminosity L 
required to have a 5a discovery in the £ ± £ ± (j/t ) sample. 



6 Final state £ + £ + £-£- 

This four-lepton signal can be produced in several decay channels of E + E~ and E ± N 
pairs, both in the case of lepton doublets and triplets, for example 

E + E~ ^£ + ZrZ, ZZ £+£-qq/is9, 

E + E~ ^ £ + Z TH/tH £~Z , Z^£ + £~,H^qq, 

E + E- -> vW + r Z/£ + ZuW~ , Z -> £ + £~ , W -»• Iv , 

E ± N -> £ ± Z£~W + , Z ^ £ + £~,W -> qq' . (35) 

For a charged singlet the signal may also result from the E + E~ decays indicated, but 
the cross section is much smaller. This signal is much cleaner than the two previous 
ones, and our event selection criteria much looser. For pre-selection we choose events 
having four charged leptons with zero total charge, two of them with transverse mo- 
mentum pt > 30 GeV and the remaining ones with > 10 GeV. For selection we 
require that opposite-charge leptons cannot be paired in such a way that both pairs 
have a mass closer to M z than 5 GeV. This latter condition is obviously included to re- 
duce the ZZnj background, but hardly affects the signals nor the rest of backgrounds. 
The number of signal and background events at the two stages of event selection is 
collected in Table [TJl The event reconstruction is performed following the procedure 
in Ref. [1], analogous to the one used for £ ± £ ± signals with missing energy: 

1. First, the two charged leptons coming from the Z boson decay are identified, 
selecting among the three possibilities the opposite sign pair £+£^ which has an 
invariant mass closest to Mz- 
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E+E- (S M ) 


37.1 


36.7 


29.8 


E+E- (NE d ) 


41.6 


41.2 


34.2 


E ± N (E M ) 


26.2 


25.8 


16.0 


E±N (NE d ) 


79.3 


78.6 


54.8 


E+E~ (E D ) 


137.6 


135.8 


111.4 


NN (NE d ) 


38.2 


38.2 


12.0 
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0.0 


E+E- (E s ) 


5.2 


5.1 


4.0 


£± iV (E D ) 


83.5 


82.6 


56.8 


NN (Z'N M ) 


62.3 


62.0 


(18.1) 


i^iV (E D ) 


0.0 


0.0 


0.0 


NN (Z'N-d) 


130.7 


128.9 


(41.5) 


tinj 


30 


30 


1 


Ztinj 


20 


20 


2 


Zbbnj 


11 


11 





ZZnj 


599 


129 


7 



Table 12: Number of £+£+£ £ events for the signals and main backgrounds with a 
luminosity of 30 fb _1 . 

2. Then, the presence of a heavy charged lepton is investigated using a plot with 
two entries per event, corresponding to the invariant mass of the Z candidate 
plus one of the remaining charged leptons, m(£ + £^£ c ) and m(£ + £^£d). 

3. Once that the location of the peak is found, we can determine which leptons 
are the E decay products by choosing between the two possibilities l+l~^l c and 
£ + £^£d, the one giving an invariant mass closest to m E . The reconstructed E 
mass m^ c is then the three-lepton invariant mass. 

We show in Fig. [8] how this procedure would work. In the upper left panel we show 
the £ + £^£ c , £t£b^d invariant mass distribution for the five sizeable signals. Those with 
a heavy charged lepton exhibit a peak, as expected, while for the rest the distributions 
are broad. On the upper right panel we show the same distribution for the background 
and also including two of the signals: the Majorana triplet (which has a resonance) and 
the Majorana singlet (which does not) to illustrate the difference. In the latter models, 
a four-lepton excess would be detected over the background but not corresponding to 
a charged resonance. Once that a peak is found, the E mass can be reconstructed, as 
shown in Fig. [8] (down, left) for the signals which display a peak. It is important to 
note that the selection between the two reconstructed mass values does not significantly 
bias the background, as seen in the lower right panel. We take the peak region as the 
interval 

280 < m^ c < 320 GeV , (36) 

and perform a kinematical cut to find the signal significances. The number of events 
at the peak are given in Table [121 For completeness, we give in parentheses the num- 
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Figure 8: Up, left: £^i^£ c , t+t^Zd invariant mass distribution for the five significant 
signals in the £ + £~£ + £~ final state (with two entries per event). Up, right: the same 
for the SM background and the background plus two of the signals. Down, left: re- 
constructed heavy charged lepton masses for several of the signals in this final state. 
Down, right: the same for the SM background and the background plus the Majorana 
triplet signal. The luminosity is 30 fb _1 . 

ber of events for the non-resonant signals. The statistical significance (neglecting the 
background uncertainty) and the luminosity required for a 5<r discovery are given in 
Tabled! 

Sq L So L 

E m 14.2 6.6 fb~ l NE d 31.3 3.0 fb~ x 
S D 52.2 1.8 ftr 1 

Table 13: Statistical significance for 30 fb _1 of the relevant signals, and luminosity L 
required to have a 5a discovery in the £ + £~£ + £~ final state. 
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7 Summary: model discrimination 



The results in the previous sections can be summarised in Table [TH which collects the 
luminosities required for 5<r discovery for each model and decay channel. The trilepton 
final state is split into the sample without Z candidates (no Z) and with a Z candidate 
(Z), and the like sign dilepton into samples without and with missing energy. The 
data in this table make apparent that the different models considered give different 
signals, and the complementarity of the channels studied make it possible to identify 
the nature of a signal eventually observed at LHC. 
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Table 14: Luminosities (in fb -1 ) required for 5a discovery for the models in the left 
column in the final states indicated. The presence of a peak in the heavy lepton pair 
invariant mass is indicated with a "P". A dash indicates an unobservable signal, or a 
discovery luminosity larger than 30 fb -1 . 

All the models except the charged lepton singlet (which will be commented below) 
give signals in the £ ± £ ± £ =F (no Z) final state, and searches in this channel are very 
sensitive to the presence of a heavy Dirac or Majorana neutrino in a singlet, doublet 
or triplet representation. For an integrated luminosity of 100 fb -1 , the approximate 
mass reach of this final state for the different models with heavy neutrinos is of 675 
GeV (S M ), 800 GeV (E D ), 850 GeV (NE d ), 850 GeV (Z'N M ) and 1 TeV (Z'N D ) for 5a 
evidence. Model discrimination is possible with the analysis of the remaining signals: 

1. Dirac / Majorana: the presence of like-sign dilepton signals without missing 
energy indicates the Majorana nature of the neutrino (models E M , Z'N M ) and 
its absence its Dirac nature (models E D , NE d or Z'Nn). 

2. Majorana triplet (E M ) / Majorana singlet (Z'N M ): these can be easily distin- 
guished because the former gives (i) four-lepton signals; (ii) much larger like-sign 
dilepton signals with missing energy; and (iii) a signal (with a peak at m E ) in 
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the £ ± £ ± £ =F (Z) final state, while the latter gives (iv) a peak in the heavy lepton 
pair invariant mass distribution, in the £ ± £ ± £ =F (no Z) and £ ± £ ± final states. 

3. Dirac singlet (Z'iV D ) / lepton doublet and Dirac triplet (NE d , £ D ): the Dirac 
singlet (i) does not give four lepton signals as the others; (ii) it has a peak in the 
heavy lepton pair invariant mass distribution. 

4. Lepton doublet (NE A ) / Dirac triplet (S D ): these the most alike signals, but 
they can nevertheless be distinguished by the presence in the case of the triplet 
of (i) a large like-sign dilepton signal with missing energy; (ii) much larger four 
lepton signals, as large as the trilepton ones; and (iii) a signal in the £ ± £ ± £ T (Z) 
final state. 

Finally, we comment on the charged singlet model. The signals in this case are unob- 
servable due to their small cross section. (Signals in the opposite-charge dilepton chan- 
nel are difficult to see even for the lepton triplet because of the large backgrounds [1,2], 
and for the charged singlet they are even smaller than in those models.) However, in a 
model with a Z' boson coupling to them, charged singlet production would be enhanced 
as it is for neutrino singlets. In such case, their presence would be characterised by: 
(i) four lepton signals with a reconstructable peak at m#; (ii) absence of £ ± £ ± £ =F (no 
Z) signals or small ones; (iii) absence of like-sign dilepton signals. 

8 Conclusions 

In order to be adequate for new physics searches, a final state signature must have 
two properties: first, it has to be sensitive to new physics, i.e. to be a possible signal 
produced in models beyond the SM; second, it must have a small background, at least 
compared to the signals expected. The trilepton final state i^i^i^ shares both features 
for heavy neutrino searches at LHC: it is sensitive to Majorana or Dirac neutrinos in 
triplet, doublet or singlet SU(2)i representations and it has a small SM background. 
The broad sensitivity of this final state, apparent with a glance at Table [HI is unique 
because other clean signals, as for example like-sign dileptons without missing energy, 
are produced by heavy Majorana neutrinos but not by Dirac neutrinos. Moreover, the 
trilepton final state is the only one in which the presence of a neutral heavy particle 
(reconstructed as a peak in the invariant mass of two opposite-charge leptons plus 
missing energy) can be established. And, with a higher luminosity, it can establish the 
triplet nature of the heavy neutrino. On the other hand, the trilepton SM background 
can be practically removed with adequate event selection criteria. For these reasons, 
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the trilepton final state can properly be considered as the golden channel for heavy 
neutrino searches at LHC. 

The like-sign dilepton final state is essential to elucidate the Majorana or Dirac 
character of a heavy neutrino. Large £ ± £ ± signals are produced in several of the models 
studied, not only from LNV processes (as in the case of a heavy Majorana neutrino) 
but also from LNC ones (as for a Dirac triplet), being the main difference the presence 
of final state neutrinos in the latter case. Therefore, we have performed two different 
analyses for the like-sign dilepton final state. In the first one our selection criteria, 
in particular the absence of significant missing energy, suppress the like-sign dilepton 
signals produced in LNC processes and efficiently isolate the true LNV ones, so that an 
eventually observed signal would correspond to a Majorana neutrino and its absence 
would indicate its Dirac character. In the second one we have required the opposite: 
large missing energy. With this novel analysis we have found that Dirac triplets give 
large and observable £ ± £ ± signals with large missing energy, while doublets and singlets 
do not. This is a notable result, because like-sign dileptons are usually regarded as a 
characteristic signature of Majorana, not Dirac fermions. (Of course, the difference in 
our analysis is the requirement of large missing energy, otherwise LNC signals would be 
suppressed.) This dilepton signal is very distinctive, and could be used to discriminate 
Dirac triplets from the other models with a heavy Dirac neutrino. 

The third final state examined, £ + £ + £~£~ , proves the presence of a heavy charged 
lepton E, with the observation of a sharp peak in a trilepton invariant mass distri- 
bution. Moreover, it can easily determine the triplet, doublet or singlet character of 
the neutrino, combined with the information from the like-sign dilepton final state. 
The charged singlet would also give a clean signal in this sample if its production is 
enhanced, e.g. by the presence of a Z' boson. The information of this four lepton 
channel is complemented by other indications, such as the observation of a peak at 
rriE in the trilepton final state with a Z candidate, or of a resonance in the heavy 
lepton pair invariant mass distribution, which help identify the underlying model with 
a luminosity not much larger than the one required for 5a discovery. Additional final 
states studied in Refs. [1,2], for example £ ± £ ± £ ± £ =F and £ ± £ ± £ ± , would provide further 
information but only at much larger integrated luminosities. 

In conclusion, in this paper we have shown that, if a positive signal of new heavy 
leptons is found, the six models considered can be discriminated already at LHC with 
the study of multi-lepton signals, without the need of heavy lepton coupling measure- 
ments at future colliders. An extra bonus of our comparative analysis is that it provides 
a guide of consistency checks to be performed (and final states to be investigated) in 
case that a multi-lepton event excess compatible with heavy lepton pair production 
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is found at LHC. We have assumed the Higgs boson to be light, as preferred by elec- 
troweak precision data, taking a mass of 115 GeV. For a heavier Higgs the results are 
expected to be similar although a dedicated analysis is required. Our study has been 
done at the level of a fast detector simulation of all the signal contributions (which is a 
non-trivial task) and the relevant SM backgrounds. The former is important, because 
there are many possible decay channels for the heavy leptons, and usually several of 
them contribute to a final state with a given charged lepton multiplicity. The latter is 
important as well. Although our simulation is not as detailed as one with a complete 
detector description, it incorporates some of the features found in real experiments, 
which parton-level analyses fail to correctly estimate. But of course, a full detector 
simulation is desirable, and it will have to be performed in order to compare the pre- 
dictions for new lepton signals with the forthcoming LHC data. 
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A Partial widths for heavy lepton decays 

We give here the partial widths for heavy lepton decays in the different models consid- 
ered, omitting the vanishing ones. For the Majorana triplet we have 

r (B + _ pw+) = (i - (i + ^ - 2^ 

327T \ m E J \ m E m 



j- j = 1 vW 

E 



, ^ g 2 .„ „m£ ( Ml\ ( Ml Mi , 
T(E + - l + Z) = -*4- \V lN 2 -f l - _f l + - 2-f = r, z 
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For the Dirac triplet, 



iw 



r{E+ -> = iw , r(iv -> rw^ + ) 

r(£?r - z-z) = r iZ , r(iv - ^z) = ^r vZ , 

r(Ef - z~#) = r lH , r(N - n H) = \v vH . 



2 



iw 



For the isodoublet the only non-vanishing ones are 

T(E + ^l + Z) = h lz , T(N^rW + ) 
T(E + -+l + H) = ±T lH , 
while for the isosinglet E we have 

r(E + - = iiw , 

T(E + ^l+Z)= l -T lZ) 
T(E+ ^l+H)= l -T lH . 
The widths for a heavy Majorana neutrino singlet are the same as for the 

T(N^l-W + ) = T lw , 

r(N^i + w-) = r lw , 

T(N — > viZ) = Y uZ , 
T(N^u l H)=T uH , 

and for a Dirac neutrino singlet they are 

T(N^rW + ) = T lw , 
T(N-^u l Z) = ^T uZ , 

T{N ^ Vl H) = \r vH . 
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